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k> , ABSTRACT 

. We report the discovery of the miUisecond pulsar PSR J2043+1711 in a search of a Fermi Large 

c3 ; Area Telescope (LAT) source with no known associations, with the Nangay Radio Telescope. The 

new pulsar, confirmed with the Green Bank Telescope, has a spin period of 2.38 ms, is relatively 
nearby (d<2 kpc), and is in a 1.48 day orbit around a low mass companion, probably a He-type 
white dwarf. Using an ephemeris based on Arecibo, Nangay, and Westerbork timing measurements, 
pulsed gamma-ray emission was detected in the data recorded by the Fermi LAT. The gamma-ray 
light curve and spectral properties are typical of other gamma-ray millisecond pulsars seen with 
Fermi. X-ray observations of the pulsar with Suzaku and the Swift /yJXJ: yielded no detection. At 
1.4 GHz we observe strong flux density variations because of interstellar diffractive scintillation, 
however a sharp peak can be observed at this frequency during bright scintillation states. At 327 
MHz the pulsar is detected with a much higher signal-to-noise ratio and its flux density is far more 
steady. However, at that frequency the Arecibo instrumentation cannot yet fully resolve the pulse 
profile. Despite that, our pulse time-of-arrival measurements have a post-fit residual rms of 2 /is. 
This and the expected stability of this system has made PSR J2043-I-1711 one of the first new 
Fermi-selected millisecond pulsars to be added to pulsar gravitational wave timing arrays. It has 
also allowed a significant measurement of relativistic delays in the times of arrival of the pulses due 
to the curvature of space-time near the companion, but not yet with enough precision to derive 
useful masses for the pulsar and the companion. Nevertheless, a mass for the pulsar between 1.7 
and 2.0 Mq can be derived if a standard millisecond pulsar formation model is assumed. In this 
article we also present a comprehensive summary of pulsar searches in Fermi LAT sources with the 
Nangay Radio Telescope to date. 

Key words: gamma-rays: general - pulsars: general - pulsars: individual (PSR J2043+1711) 
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1 INTRODUCTION 

Searches for radio pulsars coincident with Fermi Large Area 
Telescope (LAT) gamma-ray sources with no known associ- 
ations have been remarkably successful, with the discovery 
of more than 30 mil l isecond pulsars (MSPs) up to now (e.g . , 
iRansom et al. I I2OIII : iKeith et al.ll201ll : ICognard et"al]|2011l ) . 
In addition, pulsed gamma-ray emission has be en observed 
for about 20 previously-known radio MSPs (|Abdo et ahl 
I2OO9I . l2010bl : iFreire et al.ll2011al : iGuillemot et al.ll2012l ). fT 
nally, the LAT has observed gamma-ray emission from glob- 
ular clusters with properties that are co nsistent with collec- 
tive emission from populations of MSPs (|Abdo et al. I l2010dl : 
iKong et aLlboid ). MSPs are therefore an important class of 
gamma-ray sources. 

MSPs are rapidly rotating neutron stars, characterized 
by small rotational periods (P < 30 ms) and period deriva- 
tives (P<10~^^). These pulsars are thought to be "re- 
cycled", spun-up to millisecond periods by the accretion 
of matter and thus t ransfer of angular momentum from 
a binary compan ion (|Bisnovatvi-Kogan fc Kombergj 1 19741 : 
lAlpar et al.]|l982l ). More than 80% of them are in binary sys- 
tems, which makes searches for MSPs less sensitive than for 
normal pulsars. In addition, the MSPs discovered in Fermi 
LAT unassociated sources are widely distributed in Galac- 
tic latitude, whereas most radio pulsar surveys have concen- 
trated on the Galactic plane. With it s unprecedented sen - 
sitivity and localization accuracy (see lAtwood et ai]|2009l ). 
the LAT directs radio telescopes to high latitude unasso- 
ciated sources that could be unknown pulsars, missed by 
previous radio surveys. Moreover, positional uncertainties in 
the F ermi LAT Second Source Catalog (2FGL; Abdo et al.i 
|2012| ) are comparable to typical radio beam sizes, making 
radio pulsation searches very efhcient. 

Despite the numerous discoveries, 30% of the sources 
in the 2FGL catalog remain unassociated and could po- 
tentially hide unknown gamma-ray pulsars. Pulsars seen 
by the Fermi LAT are characterized by gamma-ray spec- 
tra with sharp cutoffs at a few GeV and low flux variabil- 
ity (jAbdo et a l. 2009, 2010c,). Observations of Fermi LAT 
unassociated sourc es from the First Source Catalog (IFGL 
I Abdo et al.ll2010al ) with curved spectra at the Nangay Ra- 
dio Telescope (NRT) have previously yielded the discovery 
of two radio and gam ma-ray MSPs, PSRs J2017-I-0603 and 
J2302-)-4442 ( Co gnard et all 1201 ll ). In this article we de- 
scribe a third NRT discovery, PSR J2043-|~1711, an MSP 
in a binary system located at the position of a Fermi LAT 
source with pulsar-like properties and no previously-known, 
plausible associations. 

This paper is organized as follows; in Sections 12.11 and 
12.21 we describe the radio observations of PSR J2043-I-1711. 
In Section r2.3l we discuss the timing analysis of this MSP us- 
ing the radio and the LAT data. In Sections I2.4l and l2.5l we 
describe the analysis of Fermi LAT data for the MSP, ob- 
served to emit pulsed gamma rays, and the results of Suzaku 
observations yielding no detection of X-ray emission from 
the pulsar. In Section O we discuss some of the scientific 
results, including the modeling of the observed radio and 
gamma-ray light curves of PSR J2043-I-1711 in the context 
of theoretical models of emission in the outer magnetosphere 
of pulsars, and the detection of the Shapiro delay induced by 
the companion star at superior conjunction, which allowed 



us to place constraints on the neutron star mass (Section l3.2|l 
and the detection of the proper motion and consequent lim- 
its on the distance to the pulsar (Section 13. 3p . Furthermore, 
in Section [3.31 we use the measured gamma-ray energy den- 
sity to derive tighter upper limits for the distance. Once a 
precise distance is measured, these can be used to derive a 
lower limit for the moment of inertia of the star. We present 
a final summary and discuss some scientific prospects in Sec- 
tion H 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Discovery observations 

The Fermi LAT catalog source 2FGL J2043.2-f 1711 
was already listed in the First Source Catalog as 
IFGL J2043. 2-1-1709. It has no known counterpart, and has 
spectral and variability properties that made it a plausible 
gamma-ray pulsar, with a curvature index of 12. and a vari- 
ability index of 12.2 in the IFGL catalog (see lAbdo et al.l 
l2010al . for definitions of the curvature and variability in- 
dices). A curvature index larger than 11.34 indicates that 
the source has a spectrum which departs from a pure power- 
law at the 99% confidence level, while a variability index 
larger than 23.21 implies that the source shows evidence of 
flux variability at the 99% confidence level. With its lack of 
known associations, its gamma-ray properties that are rem- 
iniscent of those of known pulsars, and high Galactic lat- 
itude (b = -15.29°), IFGL J2043.2-hl709 satisfied aU the 
criteria used for selecting the sources eligible for the orig- 
inal search for pulsars in Fermi LAT unassociated sources 
at Nangay that l ed to the discovery o f PSRs J2017-I-0603 
and J2302-I-4442 (|Cognard et aLlbOllI ). except for the spa- 
tial localization: the semi-major axis of the IFGL source 95% 
confidence ellipse (S95) was larger than the conservative cut 
of 3' used for the latter search in order for sources to be well 
covered by the Nangay beam, which has a width at half max- 
imum of 4' in right ascension and 22' in declination. With 
additional data, the source localization improved (S95 = 2.9' 
in the 2FGL catalog), making 2FGL J2043. 2-1-1711 an ex- 
cellent candidate source for radio pulsar searches with the 
NRlQ. 

A first 1-hr observation of the source was made 
on 2009 November 21 at the NRT, using the mod- 
ified Berkeley-Orlcans-Nangay (BON) instrumentation 
(Cognard & Thourcau 2006) at 1.4 GHz, with a 512 x 0.25 
MHz incoherent filter bank sampled every 32 /is. The data 
were dedispersed in ~2000 dispersion measure (DM) values, 
up to 1244 pc cm""^, and processed using acceleration and 
single pu lse search techniques as provided by the Presto 
package (R ansom e t al. 2002). No pulsations were observed 
in this initial observation. Similarly, the analysis of a second 
observation taken on 2009 December 2 yielded no detection. 
Nevertheless, a 19cr candidate with a rotational period of 
2.379 ms and a DM of 20.7 pc cm"'' was observed in the 
third observation, taken on 2009 December 12. Seven sub- 
sequent observations of the Fermi LAT source at 1.4 GHz 

^ In addition to the successful search strategy just described, a 
number of other NRT observations of other Fermi LAT sources 
are documented in the Appendix. 
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with the NRT yielded no re-detection, casting doubts on the 
presence of a pulsar in this Fermi LAT source. The 2.379 
ms candidate was finally confirmed with observations made 
at the GBT telescope at 350 Mffz with the GUPPI back- 
enjf during a survey o f Fermi LAT unassociated sources 
(see 
327 



Hessels et al.ll20l"lh . and at the Arecibo telescope at 
MHz with the Wide-b and Arecibo Pulsar Processors 
(WAPPsi lOowd et al.ll200Cl ). Substantial accelerations of the 
rotational period across the confirmation observations were 
measured, indicating that the pulsar is in a binary system. 



2.2 Arecibo Observations 

Radio light curves of PSR J2043-I-1711 recorded with the 
Arecibo telescope are shown in Figure [T] At 1.4 GHz the 
pulsar is observed to exhibit dramatic radio fiux variations, 
explaining the several unsuccessful attempts to confirm the 
pulsar at Nangay following the discovery. Even with Arecibo 
it can be a difficult object: during one 30-min observation 
the pulsar was not visible in one of the WAPPs, centered at 
1.46 GHz and with 50 MHz of bandwidth. Assuming a sys- 
tem temperature of Tsys = 25 K, a gain of G = 10 K Jy~^, 
a pulse width of 6.25% (8 bins above average out of a total 
of 128), and a signal-to-noise ratio (SNR) smaller than 3, 
this corresponds to an instantaneous upper limit on the ra- 
dio flux density of 5.1 /iJy. However, on most occasions the 
pulsar is detectable with Arecibo and the high SNR light 
curve was recorded during a bright scintillation state. The 
1.4-GHz light curve is complex, with several pulsed com- 
ponents. A total of 4 useful observations of the pulsar at 
1.4 GHz have been carried out with Arecibo, between 2010 
November 20 and 2011 August 24. The average radio fiux 
density for these observations was of the order of 10 ^Jy. 

At 327 MHz the flux density is far more steady, and for 
that reason we carried out the bulk of the Arecibo timing at 
this frequency. The pulsar was observed 64 times between 
2010 July 17 and 2011 August 25 with average integrations 
of 35 minutes using the 327 MHz receiver (G = 11 K Jy~^, 
Tsys = 113 K). For most observations we use the 4 WAPP 
spectrometers in parallel. Each WAPP makes a 3-level dig- 
itization of the analog voltages over a 12.5 MHz band for 
both linear polarizations, autocorrelating these for a total 
of 512 lags. The data are then integrated for ts = 64 ^s 
and the orthogonal polarizations added in quadrature are 
written to disk. Their bands are centered at 308.25, 320.75, 
333.25 and 345.75 MHz, and together they cover the fuU 50 
MHz band provided by the receiver. For all Arecibo observa- 
tions, the lags were Fourier transformed to generate power 
spectra. These were dedispersed at the nominal DM of the 
pulsar and folded modulo its spin period using the Presto 
pulsar software package, generating pulse profiles that are 
then stored for later analysis. Ten Arecibo observations at 
345 MHz were averaged to derive a radio fiux density of 
'S'345 ~ 1-2 mjy with a standard deviation of 0.2 mjy, as- 
suming a pulse width of 6.25%. 

Note that this pulse width is overestimated because of 
instrumental limitations: with the WAPPs the pulses are 
seen with an effective width at half maximum dt given by 



https://safe.nrao.edu/wiki/bin/view/CICADA/NGNPP 
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Figure 1. Multi-wavelength light curves of PSR J2043+1711. 
The bottom panel shows radio profiles recorded at the Arecibo 
telescope at 345 MHz (40-min integration), at the Westerbork 
Synthesis Radio Telescope at 345 MHz (3.8-hr integration) and 
at the Arecibo telescope at 1.4 GHz (1-hr integration). The third 
panel shows a 40-bin gamma-ray profile obtained by selecting 
events recorded by the Fermi LAT within 5° from the pulsar and 
with energies above 0.1 GeV, and weighting each event by its 
probability of originating from the pulsar. The fitting functions 
are shown as dashed lines in the first cycle. The top two panels 
show non-weighted Fermi LAT light curves for events recorded 
within 0.8° of the pulsar, with energies between 0.1 and 1 GeV, 
and above 1 GeV, respectively. No gamma-ray features signif- 
icantly narrower than the bin width used here were observed. 
Horizontal dashed lines indicate gamma-ray background levels. 
Two rotations are shown for clarity. 



the sum in quadrature of the intrinsic width at half maxi- 
mum dti , the sampling time dts and the dispersive smearing 
dtd- At an observing frequency v =0.345 GHz, we have 
dti ~ 75 — 80 /is (see below). The sampling time is 64 /is, 
and with a bandwidth of 12.5 MHz for each WAPP, and 
512 spectral channels across that band width, the dispersive 
smearing dtd (see e.g. Equation (5.2) of iLorimer fc Kramer! 
I2OO5) is calculated to be ^ 120 /is. Therefore, we have an 
effective width of ~ 156 /is, twice as much as the intrinsic 
pulse width. The int rinsic pulse profile was dete rmined using 
the PuMa2 backend (jKaruppusamv et al.ll2008l ') at the West- 
erbork Synthesis Radio Telescope (WSRT), in the Nether- 
lands, which is capable of coherent dedispersion for a total 
bandwidth of 80 MHz centered at 345 MHz (see Figure [TJ . 
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2.3 Timing analysis 

The best detections of the pulsar at each instrument and 
frequency were used to derive "standard" pulse profiles. For 
the dominant Arecibo 327 MHz dataset, we made 1 time of 
arrival (TOA) per WAPP for every 500 s of observations on 
average, by cross-correlating the puls e profiles with the stan- 
dard profiles in the Fourier domain (jTavloil [l992h . We also 
extracted 32 TOAs from the 1.4 GHz Arecibo observations. 
This resulted in a total of 1029 Arecibo TOAs. In addi- 
tion, 18 NRT TOAs have been recorded at 1.4 and 1.6 GHz 
from 2010 August 16 to 2011 Aug ust 15 using the proce dure 
and instrumentation described in lCognard et al.l(|201ll ). Six 
TOAs were recorded at 345 MHz with the WSRT between 
2011 March 18 and 2011 August 20. Finally, gamma-ray pul- 
sations from PSR J2043-I-1711 were detected in the Ferrm 
LAT data after the first months of radio timing following 
the discovery (see Section [2]4]), which allowed us to measure 
T OAs using the m aximum likelihood techniques described 
m iRav et al.l (|201ll ) and recover phase-coherence across the 
entire Fermi LAT dataset. A total of 13 gamma-ray TOAs 
with average uncertainty 25.3 /is and corresponding to at 
least 3(T detections were extracted between 2008 Septem- 
ber 12 and 2011 June 10 by selecting photons with ener- 
gies greater than 0.5 GeV and with reconstructed directions 
found within 1° of the pulsar. 

We carried out subse quent TOA analy ses using the 
TEMP02 software package (|Hobbs et aLlbOOq ). For the con- 
version of Terrestrial Time (TT) TOAs to Coordinated 
Barycentric Time (TCB) we used th e DE/LE 421 solar 
system ephemeris (|Folkner et al.l [20091 ). The differences be- 
tween observed and predicted barycentric TOAs (the tim- 
ing residuals) were weighted in the fit according to the es- 
timated uncertainty of each TOA. The resulting timing pa- 
rameters and their la uncertainties are presented in Table 
[U The orbit of PSR J2043+1711 has very lo w eccentricity 
therefore we us ed the "ELLIH" orbital model (|Freire fc Wexl 
I2OIOI : llTange et al . 2001) to model it. This yields Keplerian 
(semi-major axis of the pulsar orbit projected along the 
line of sight, x, orbital period Pj,, epoch of ascending node 
To, eccentricity e and longitude of periastron, u)) and post- 
Keplerian (orthometric amplitude and orthometric ratio 
?) parameters that are weakly correlated with each other. 

This ephemeris describes the TOAs well, with a reduced 
of 2.1 for 1044 degrees of freedom. In order to estimate the 
timing parameters with realistic uncertainties, we adjusted 
the uncertainty estimates of each individual timing dataset 
using EFAC parameters so that x^/'^frcc is equal to 1 in ev- 
ery case. The uncertainties quoted in Table[l]were derived in 
this fit by TEMP02, except where stated otherwise. The 1029 
Arecibo TOAs have a post-fit rms uncertainty of 2.13 /^s, 
despite the aforementioned smearing caused by the instru- 
mentation used to date. For this reason PSR J2043+1711 
has been added to the Int ernational Pulsar Timing Array 
(IPTA) (|Hobbs et aLlbOld l. 

There is clearly scope for further improvement in the 
timing precision of this object, given the fact that the pulse 
width measured with the WAPPs is twice as wide as the 
intrinsic pulse width. This implies that, with a broadband 
coherent dedispersion instrument working at Arecibo, the 
measured peak flux density should be about twice as large 



and the rise time about half as long, implying about four 
times the current timing precision. 

Note that the DM value and its uncertainty were mea- 
sured independently from the analysis described here. To 
measure the DM we built a dataset of 308.25, 320.75, 333.25 
and 345.75 MHz Arecibo TOAs by cross-correlating the indi- 
vidual pulse profiles with a single standard profile, to prevent 
any phase shifts caused by the usage of different template 
profiles. This dataset was then fitted for the DM, yielding 
the best-fit value and corresponding uncertainty listed in 
Table □ 

Because the radio timing data cover only slightly more 
than a year, no proper motion could be measured with these 
data alone. However, with the gamma-ray TOAs covering 
approximately three years, we could measure a significant 
proper motion (see Table [l]). It is clear that the proper 
motion measurement depends strongly on the gamma-ray 
timing. We therefore checked the proper motion values and 
their associated uncertainties by studying the influence of 
the LAT timing on the measurement. We flrst repeated the 
analysis described above using different numbers of gamma- 
ray TOAs, from 5 to 25 with a step of 1. This yielded average 
proper motion values that are consistent with those listed in 
Table [1] within standard deviations on cos ((5) and ^is of 
0.1 and 0.4 mas yr~^, respectively. Also, we made 1000 re- 
alizations of a Monte-Carlo simulation in which gamma-ray 
TOAs were generated based on the timing solution given in 
Table [U and were then perturbed so that the residuals have 
the same rms as those of the actual gamma-ray TOAs (~ 
30 /.ts). The uncertainties were flnally shuffled from the ac- 
tual gamma-ray TOAs. Again, this yielded average /Iq cos((5) 
and /ia values that are compatible with the values reported 
here, within standard deviations of 0.6 and 0.7 mas yr~^, 
respectively. Combining the standard deviations obtained 
from these two studies, we estimate that the gamma-ray 
TOAs introduce systematic uncertainties on /i^ cos (5) and 
of approximately 1 mas yr~^. 



2.4 Gamma-ray analysis 

The analysis of the LAT data was done using the Fermi 
Science Tool^ (STs) v9r23pl. We selected events recorded 
between 2008 August 4 and 2011 July 21, with reconstructed 
directions within a 20° x 20° region centered on the pulsar 
position, energies between 0.1 and 100 GeV, and zenith an- 
gles below 100°. We further selected "Source" class events 
of the P7_V6 instrument response functions (IRFs), and re- 
jected times when the rocking angle of the telescope ex- 
ceeded 52° and when the Earth's limb infringed upon the 
region of interest (ROI) . The gamma-ray events were phase- 
folded using the Fermi plug-in distributed with Tempo2 
(|Rav et al.ll201ll ) and the ephemeris given in Table [T] 

To measure the spectral properties of the pulsar we 
fltted sources in the ROI using a binned maximum likeli- 
hood method, using the pyLikelihood module provided with 
the STs. The spectral parameters of the 54 2FGL catalog 
sources within 20° of the pulsar were included in the model. 
PSR J2043+1711 and the eight other pulsars in the held of 



^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 
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Table 1. Measured and derived parameters for PSR J2043+1711. Numbers in parentheses are the nominal Icr TEMP02 uncertainties in 
theleast-significant digits quoted. The distance was estimated using the NE2001 model of Galactic free electron density (ICordes &: Laziol 
I2OO2I) . We assumed an uncertainty on this distance estimate of 20%. Using this distance estimate and the measured proper motion, 
/iy, we calculated the period derivative corrected for the Shklovskii effect (Shklovskii IQTOj), -Pcorn and used that value to derive E, 
and Bi^c- Note that E, Bs, Bj^q and 7? were calculated assuming a moment of inertia / of 10'*^ g cm^. For proper motion parameters 
and for the total proper motion, the first quoted uncertainties are the la statistical uncertainties from TEMP02 and the second ones are 
systematic (see Section 12.31 for details on the calculation of the systematic uncertainties) . For gamma-ray parameters the first quoted 
uncertainties are statistical and the second are systematic, and correspond to the differences observed when doing the spectral analyses 
with the P6_V3 IRFs and associated diffuse models. 



R.A. (J2000) 20^43^20=88309(5) 

Decl. (J2000) +17° ll'28y 948(1) 

Rotational period, P (ms) 2.37987896026(4) 

Apparent period derivative, P (10~-^^) 5.24(2) 

Proper motion in right ascension, cos((5) (mas yr~^) . . . —7 ± 1 ± 1 

Proper motion in declination, fig (mas yr~^) —11 ± 2 ± 1 

Epoch of ephemeris, Tq (MJD) 55400.00019 

Dispersion measure, DM (cm-3 pc) 20.70987(3) 

Orbital period, Pj, (d) 1.482290809(2) 

Projected semi-major axis, x (It s) 1.6239614(1) 

Epoch of ascending node, Tasc (MJD) 55253.8038503(6) 

rj = esinoj (10"^) -2.1(1) 

K = ecosu (10-^) -2.6(1) 

ha 0.63(7) 

T (/is) 0.87(5) 

Span of timing data (MJD) 54729.1 — 55798.2 

Number of TOAs 1066 

RMS of TOA residuals (ps) 2.1 

Solar system ephemeris model DE421 

Units TCB 

Flux density at 327 MHz, ^327 (mjy) 1.2 ± 0.2 



Derived parameters 

Orbital eccentricity, e (10~®) 

Mass function, / {Mq) 

Minimum companion mass, nic (Mq) 

Galactic longitude, Z (°) 

Galactic latitude, (°) 

Distance inferred from the NE2001 model, d (kpc) 

Total proper motion, fiT (mas yr" ^) 

Shklovskii-corrected period derivative, Pcorr (10~^^) 

Spin-down luminosity, E (10'^'* erg s" ^) 

Surface magnetic field strength, Bs (10^ G) 

Magnetic field strength at the light cylinder, B^c (10* G) 

Light curve parameters 

First peak position, $1 

First peak full width at half-maximum, FWHMi . . 

Second peak position, <I>2 

Second peak full width at half-maximum, FWHM2 

Radio-to-gamma-ray lag, S 

Gamma-ray peak separation, A 

Gamma-ray spectral parameters 

Spectral index, F 1.4 ± 0.1 ± 0.4 

Cutoff energy, Ec (GeV) 3.2 ± 0.6 ± 1.0 

Photon flux (> 0.1 GeV) (10"* cm-^ s"!) 2.8 ± 0.3 ± 0.9 

Energy flux (> 0.1 GeV), F^, (10"" erg cm'^ s"!) 2.8 ± 0.2 ± 0.3 

Luminosity, L^ / fn (10^* erg s^^) (1.0 ± 0.1 ± 0.1) X (cZ/1.8 kpc)^ 

Efficiency, J? / fn (1.0 ± 0.1 ± 0.1) X (d/1.8 kpc)^ 



3.4(1) 
0.00209287(8) 
> 0.173 
61.92 
-15.31 
1.8 ±0.4 
13 ± 2 ± 1 
3.6 ±0.5 
1.1 ±0.2 
9.3 ±0.7 
6.3 ±0.5 



0.17 ±0.01 
0.10 ±0.02 
0.58 ±0.01 
0.10 ±0.03 
0.18 ±0.01 
0.41 ±0.02 
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view were modeled as exponentially-cutoff power laws of the 
form: 




In this equation, TVo denotes a normalization factor, 
r is the photon index, and Ec is the cutoff energy of 
the spectrum. The extragalactic diffuse emission and the 
residual instrument background were modeled jointly us- 
ing the iso_p7v6source template, while the Galactic diffuse 
emission was modeled using the gaL2yearp7v6_vO mapcube. 
The spectral parameters of the 12 sources within 10° of 
PSR J2043-I-1711 were left free in the fit, as were the nor- 
malization factors for the diffuse components. The nearest 
source in the catalog, 2FGL J2031. 0-1-1938, is located ~ 3.7° 
away and its gamma-ray energy flux above 0.1 GeV is 3 
times smaller than that of PSR J2043-I-1711. The measure- 
ment of spectral properties for the latter object is therefore 
weakly affected by the neighboring sources. The gamma- 
ray spectral energy distribution of PSR J2043-f 1711 for an 
exponentially cutoff power law is shown in Figure [2] The 
best-fit spectral parameters, as well as the integrated pho- 
ton and energy fluxes above 0.1 GeV, are listed in Table [T] 
With the measured energy flux F^, we derived the gamma- 
ray luminosity Lj = ATvfnFjd^ and efficiency of conversion 
of spin-down power into gamma-ray emission, ri = L-y/E, as- 
sumin g a geometrical correction factor fn fsee lWatters et al] 
(|2009( ) for the definition) of 1, which is typical under 
outer- magnetosph eric gamma-ray emission models (see e.g. 
I Venter et aLll2009^ . The best-fit spectral parameters for the 
sources within 10° were compatible within statistical and 
systematic uncertainties with the values published in the 
2FGL catalog, with the exception of 2FGL J2035.4+1058, 
located 6.5° away from PSR J2043-I-1711, for which our flux 
estimate is larger than the 2FGL flux by four standard devi- 
ations. The 2FGL source J2035. 4-1-1058, which is associated 
with the blazar PKS 2032-1-107, is nevertheless flagged in 
the catalog as being variable. Since the time interval consid- 
ered in our analysis was longer than in the 2FGL catalog, it 
is not surprising that the flux we measured differs from the 
2FGL one. We checked the best-flt spectral parameters for 
PSR J2043-I-1711 b y fltting the data with the pomthke hke- 
lihood analysis tool (Kerr 2011b), and found results that are 
consistent with those listed in Table [T] within uncertainties. 

Using the full spectral model obtained with the like- 
lihood analysis and the tool gtsrcprob, we calculated the 
probabilities that each gamma-ray event originates from the 
pulsar. The probability-weighted gamma-ray light curve of 
PSR J2043-I-1711 above 0.1 GeV and for events found within 
5° of the p ulsar is show n in Figure [T] The weighted _ff-test 
parameter (|Kerij[2011al ) is 433.5, corresponding to a pulsa- 
tion signiflcance of 18.5cr. The upper two phase histograms 
in Figure [T] show gamma-ray light curves for events found 
within 0.8° of the pulsar position, in different energy bands. 
The background levels in these light curves have been calcu- 
lated by summing the probabilities that events have not been 
pro duced by the MSP, as de scribed in .Cognard et al.l (1201 ll ) 
and iGuillemot et al.1 l|2012h . We verifled that the weighted 
light curve did not show indications for emission features 
significantly narrower than the bin width used, by analyz- 
ing the same light curve with five and ten times the number 




Energy (GeV) 

Figure 2. Phase-averaged gamma-ray energy spectrum for 
PSR J2043-fl711. The black hne shows the best-fit model ob- 
tained by fitting the Fermi LAT data with a simple exponentially 
cutoff power-law form (see Section 12.41 for details), while dashed 
error lines indicate la errors. Data points are derived from like- 
lihood fits of individual, variable- width energy bands defined by 
the requi rement that the pul sar be detected with a Test Statistic 
(TS; see lMattox et aLllToQ^) of at least 50. In these bands the 
pulsar is modeled with a simple power-law form. An upper limit 
was calculated for the last energy band as the pulsar was not 
detected with enough significance in that band. 



of bins, finding no statistically significant component other 
than those seen in Figure [1] 

We measured the phase and full width at half max- 
imum FWHMi of each gamma-ray component by fitting 
the two peaks in the weighted light curve with Lorentzian 
functions. The best-fit parameters are listed in Table [TJ 
along with the radio-to-gamma-ray lag, 5 = $i — where 
= 0.99 is the phase of the maximum of the 345 MHz ra- 
dio profile, and the separation between the two gamma-ray 
peaks, A = $2 — $1. 

The gamma-ray light curve shape of PSR J2043-I-1711 
is similar to those of other gamma -ray MSPs (see e.g. 
lAbdo et al.ll2009l : ICognard et al.ll201ll '). suggesting that the 
gamma-ray emiss ion is produced at high altitudes in their 
magnetospheres l|Venter et al.l |2009| ). Likewise, the spec- 
tral properties of t he MSP are simil ar to those of other 
gamma-ray pulsars (|Abdo et al.ll20103 ). The large efficiency 
rj ^ 100% indicates that the distance is very likely over- 
estimated by the NE2001 model. In Section [3] we however 
discuss the possibility to use this large efficiency value to 
constrain the moment of inertia of the star, if the actual 
pulsar distance is close to the DM distance. 



2.5 X-ray analysis 

On 2010 May 3, the Fermi LAT unassociated source 
IFGL J2043. 2-1-1709 was observed by the Suzaku X-ray ob- 
servatory for 18 ks, as part of an unassociated sources ob- 
servation campaign. A description of this campaign and the 
data reduction, as well as an X-ray image of the field of 
view aroun d IFGL J2043. 2-1 - 1709 as seen with Suzaku can 
be found in iTakahashi et al.l (|2012l '). 

No significant X-ray source can be seen at the po- 
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sition of PSR J2043+1711. Assuming a power- law model 
with a photon index of 2, the X-ray upper limit for 
PSR J2043-I-1711 between 0.5 and 8 keV is calculated to 
be - 4.7 X 10"^* erg cm'^ s"^ (90% confidence). 

We also checked the available Swift /XI{T observations 
of PSR J2043-I-1711, and found two, totaling - 10.5 ks of 
data. No X-ray sources were detected within 2' of the pulsar 
position, thus confirming the Suzaku results. 



3 DISCUSSION 

3.1 Light Curve Modeling 

We have fit the radio and gamma-ray light curves of 
PSR J2043-I-1711 to geometric simulations (assu ming the 
vacuu m retarded dipole magnetic field geometry of lDeutschI 
Il955h using a maximum likelihood technique. The gamma- 
ray light curve s were fit with the two-pole c austic (TPC; 
Dvks fc RudaS |2003| ) and outer-gap (OG; ICheng et all 
19861 ) models. For our purposes the TPC model is taken 
to be a geometric realiz ation of the slot-gap model 



l|Muslimov fc Hardindl2004h . The radio light curves were fit 
with a hollow-cone be am, core beam, an d cone plus core 
beam models following IStorv et al.l (|2007l ). 

The gamma-ray light curve (30 bins per rotation) used 
for the fits was constructed using all events from the data 
described in Section r2.4l with reconstructed directions within 
0?8 of the pulsar radio position. Additionally, we used either 
the 1.4 GHz or 345 MHz radio profile (also in 30 bins) but 
report only results from fits with the former as the latter 
profile is known to be too wide; however, fits using the 345- 
MHz profile were used in estimating systematic biases in our 
procedure as discussed below. 

We have used simul ations with th e same parameters 
(P — 2.5 ms) as those in lCognard et all (|201ll ) except that 
we have a r esolution of 2.5% of the polar cap opening angle 
(Opc ~ ^J2-kRn s / Pc) in gap width. We scanned over the 
parameter phase space in order to find the best-fit for each 
model as given in Table [2l 

To account for the fact that these models are rel- 
atively simple and the likelihood surfaces can be steep 
around maxima, implying unrealistically small uncer- 
tainties, we have rescaled the likelihood differences by 
rifree/ [2 * (— \n(Cmax))\, where Cmax IS the maximum likeli- 
hood value and nfroo is the degrees of freedom for the given 
fit. Assuming the log-likelihood differences follow a distri- 
bution this results in the best fit corresponding to a reduced 
= 1. The best-fit uncertainties given in Table [2] are 68% 
confidence level. The two widths reported are for the accel- 
erating {wacc) and emitting (wem) gaps, these are the same 
for the TPC model. 

Our models are relatively simple and, in the case of the 
radio profiles, do not contain as many components as are im- 
plied by the data. Additionally, the geometry used for the 
magnetic field cannot be correct as the magnetosphere will, 
to so me extent, be filled with charges (|Goldreich fc JulianI 
1 19691 ) which will distort the gap geometry. Therefore, we 
have attempted to estimate systematic biases in the best-fit 
geometries reported in Tabled To do this we refit the data 
while varying the radio uncertainty by a factor of 2, varying 
the gamma-ray background estimate by 5%, using the 1.4 
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Figure 3. Data and best-fit Hght curves using the gamma-ray 
data described in Section [STT] (top) and the 1.4 GHz radio profile 
(bottom). Models corresponding to the TPC fit are shown as 
dash-dotted lines (solid pink lines online) and for the OG fit as 
dashed lines (solid green lines online). 



GHz radio light curve in 60 bins, and performing fits with 
gamma-ray light curves corresponding to the energy ranges 

0. 1-1 GeV and > 1 GeV. None of these changes strongly af- 
fected the gap width parameters. Changes in geometry were 
typically ^8°; however, the OG best-fit geometry changed 
by 25 - 40° in one parameter for both the core and cone 
beam fits when the radio uncertainty was doubled and when 
fitting only the 0.1-1 GeV gamma-ray light curve. 

The data and best-fit model light curves are shown in 
Figure 13] using the 1.4 GHz radio profile and the hollow-cone 
beam model. While the likelihood formally favors the TPC 
model, both gamma-ray models reproduce the qualitative 
features of the observed gamma-ray light curve (gamma-ray 
peak separation and bridge emission) well, and the radio 
model is not optimal. The OG model therefore cannot be 
ruled out. Fits with both the TPC and OG model do find 
solutions with C, near 80° which is to be expected if the spin 
and orbital axes have (nearly) aligned over time. 

It should be noted that the likelihood favors the core 
beam model in all cases. Nevertheless, under these fits we 
obtain impact parameters (/? = C, — a) oi order 20° which 
would imply a very faint radio fiux for a beam falling ofi^ as a 
Gaussian away from the magnetic axis. Such solutions thus 
seem less likely. Note, however, that no polarimetric data 
exist for this MSP to conclusively confirm or rule out the 
presence of a core beam. 

From our models we have estimated /n (described in 
Section 2.4) for each model and provide estimated uncertain- 
ties (see Table [2ll. In all cases the predicted /n is less than 

1, although relatively close. This leads to corrected gamma- 
ray efficiencies of = (0.8 ± 0.1 ± 0.1) x (ci/1.8 kpc)^ for the 
TPC model and = (0.9l5;;3l2;3) x (d/1.8 kpc)^ for the OG 
model. 



3.2 Component masses 

From the projected semi-major axis x and the orbital 
period Pb we calculate the mass function f{mp,mc) = 
{mcsmif/{mp + mcf = {4n^e'x^)/{GMQPi) ~ 2.1 x 
10~^Mq, where rup, rric and i are the pulsar mass, the com- 
panion mass and the orbital inclination. Assuming a pul- 
sar mass rrip of 1.4 Mq, we derive a minimum companion 



8 L. Guillemot et al. 



Table 2. Best-fit parameters from fight curve modeUng (see Section l3.ip . For the gap width parameters the maximum size is 0.10. Values 
of 0.0 are unphysical and should be taken to mean that the best-fit width is less than our resolution of 0.025. For the OG models the 
width of the emitting gap is constrained to be no more than one-half the accelerating gap size. 
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mass of rric > 0.173 Mq. The mass of the companion could 
be significantly higher for much lower orbital inclinations. 
However, that possibility can be excluded by our detection 
of relativistic time delays in the TOAs caused by space- 
time curvature in the vicinity of the companion star, com- 
monly known as the Shapiro delay (Shaoiro 1964j); these are 
displayed in Figure |4] This detection has high significance: 
the orthometric amplitude (^3) measurement is nine times 
larger than the la uncertainty. This detection implies that 
the orbital inclination must be high. The constraints on rric, 
rUp and sin i introduced by our detection of the Shapiro de- 
lay are depicted in Figure [S] The mass function and range 
of companion mass values suggest that the companion star 
is likely to be a He- type white dwarf (WD). 

The second Shapiro delay parameter, the orthomet- 
ric ratio (?), is not yet precise enough to determine as- 
trophysically meaningful values for nic, mp and sini. We 
can, ne vertheless, make an estim ate of rUp if we assume 
that the iTauris fc Savoniid (|l999l ) relation between ntc and 
Pb applies, as it does for all known MSP/He- type- WD 
syste ms with precisely measu red masses (see e.g. Figure 
2 of Ivan Kerkwiik etaP I2OO5I I. For the orbital period of 
PSR J2043-I-1711, we derive 0.20 < rUc < 0.22 Mq, which 
is very close to the best value derived from the Shapiro delay 
measurements as can be seen from Figure [S] In this case, 
the existing constraints imply 1.7 < rUp < 2.0 Mq and 
i — 81.3 ± 1.0 degrees (that is, a nearly edge-on configu- 
ration). Thus the mass of this pulsar ap pears to be located 
between the masses of P SR J1903-I-0327 (iFreire et al.l2011bl ') 
and PSR J1614-2230 (iDemorest et al.ll2010l ). which define 
the high end of the millisecond pulsar mass distribution. 



3.3 Constraints on the distance and the moment 
of inertia 

With the current set of TOAs we were able to measure the 
pulsar's proper motion, finding fia cos(5) = — 7± 1 mas yr"^ 
and fj,s = —11 ±2 mas yr"^, corresponding to a total proper 
motion of /^t = 13 ± 2 mas yr"^. For a distance of 1.8 kpc, 
this gives a transverse velocity Vr ~ 110 km s"^, a value 
that is relativel y typical among G alactic disk millisecond 
pulsars (see e.g. lHobbs et al.ll2005l ). 

This transverse motion induces a const antly changing 
Doppler shift first noted bv lShklovskiH (|l970l '). which makes 
the apparent P value greater than the intrinsic one by Ps ~ 
2.43 X 10~2i s~i Pdfi^, where P is the rotational period 
and d is the distance. The NE2001 model of Galactic free 
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Figure 4. Arrival time residuals for the 0.3 GHz Arecibo TOAs 
of PSR J2043-I-1711, as a function of orbital phase. The timing 
residuals were binned in phase, with 20 bins per orbit. TOA un- 
certainties were taken into account when binning the timing resid- 
uals. Bottom: residuals obtained with the full timing model listed 
in Tablc[T] Top: residuals for best-fit orbital model that does not 
take into account the Shapiro delay. The dashed line shows the 
theoretical prediction for the detectable part of the Shapiro delay 
(which is not absorbed by the fitting of the Kepleria n parameters 
X and e) given by Eq. (19) in iFreire fc Wexl hoid) and the hs 
and ? parameters in Table [l] 



electron density distrib ution places the pulsar at 1.8±0.4 kpc 
(ICordes fc Lazidl2002l ). For this distance and proper motion. 
Pa is found to represent ~30% of the measured (apparent) P 
value. The Shklovskii-corrected period derivative Pcorr and 
the derived pulsar properties are listed in Table [1] 

An upper limit on the distance can be derived by assum- 
ing that the Shklovskii effect accounts for all of the apparent 
P and that the currently measured proper motion is correct. 
Doing so, we find dmax = 5.7 kpc (see Figure [6]). It is clear, 
however, that the pulsar must be at a distance from the Solar 
System that is significantly smaller than dmax- The reason 
for this has already been briefiy mentioned in Section 12.41 
namely the very high implied gamma-ray efficiency: for dis- 
tances larger than 1.8 kpc, an efficiency larger than 100% 
is required to produce the gamma-ray flux detected by the 
LAT. 

We can derive more realistic upper limits for d based on 
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Figure 5. Constraints on some physical parameters (companion mass rUc, pulsar mass rrip and orbital inclination i) of the 
PSR J2043+1711 binary system. The purple curves enclose regions consistent with the nominal values and la uncertainties of the 
Shapiro delay parameters (solid) and if (dashed). Left: nic-cosi plot. The gray region is excluded by the condition nip > 0. The 
black solid curve is a contour level of the 2D probability density function (PDF) that encloses 68.3% of the total probability. Bight: 
rric-mp plot. The gray region is excluded by the condition sini < 1. The black solid curve encloses 68.3% of the total probability in this 
region of the plane, and is not a translation of the black contour curve in the left-hand plot. Top, right marginal plots: the solid black 
lines show the ID probability distribution functions for cosi, and m,c. For a l l plot s, the red contours represent sections of the same 
2D PDF for which 0.20 < nic < 0.22 AIq according to the lTauris fc Savoniij ll 19991 ) model based on the orbital period measured for 
PSR J2043+1711. 



our Fermi LAT measurement of the energy density G and 
imposing the condition 77 < 1. A comparison of the expres- 
sions for E and then yields the following inequality: 



/ > 



(2) 



For the current assumptions of / = 10*^ g cm'^ and /n — 
1, this inequality only holds for distances smaller than 
1.8 kpc (Figure [6|. For the best-fit fn factors obtained 
from the modeling of radio and gamma-ray hght curves of 
PSR J2043-hl711 with the TPC and OG models and with 
the hollow-cone beam radio model, the inequality holds for 
d < 1.92 and 1.83 kpc, respectively. The larger distances 
within these ranges are comparable to the distance predicted 
by the NE2001 model. 

Improved radio timing of this pulsar might eventually 
allow a precise measurement of the timing parallax, and 
therefore of the distance d, which according to Eq. [2] pro- 
vides a direct lower limit for I. Again, such limits are dis- 
played in Figure [S] As was mentioned in Section 12.41 the 
NE2001 model likely overestimates the distance. However, 
it is clear that if the model is correct or the pulsar is farther 
away, then its moment of inertia must be large. A high lower 
limit on the moment of inertia would help co nstrain equa- 
tions of state (EoSs) for super-dense matter (jWorlev et al.l 



|2008| ). especially if combined with a precise measurement 
for the mass of the pulsar, another likely consequence of the 
improved timing. 



4 CONCLUSION AND PROSPECTS 

We have reported the discovery of an MSP with the Nangay 
Radio Telescope at the position of an unassociated Fermi 
source, PSR J2043-I-1711. The pulsar is the third MSP to 
be discovered at Nangay in a Fermi source, aft er PSRs 
J2017+0603 and J23024-4442 (|Cognard et al.ll201lh . The ra- 
dio pulsar is found to be responsible for the gamma-ray emis- 
sion observed by Fermi, and its properties (rotational pe- 
riod, spin-down luminosity, distance, gamma-ray light curve 
and gamma-ray spectrum) are relatively common among 
known gamma-ray MSPs. 

Of the pulsars discovered in Ferm i unassociated sources 
that have been published to date ( Cognard et all I2OIII : 
iRansom et all 1201 ll : iKeith et al.ll201ll ). PSRs J2043-H711 
and J2017-I-0603 are the two systems with the best timing 
precision. This happens because they have very sharp fea- 
tures in their pulse profiles. Both objects are relatively faint, 
but they are well within the region of the sky detectable 
by the Arecibo 305-m telescope, which greatly compensates 
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Figure 6. Upper limits on the distance as a function of the mo- 
ment of inertia, imposed by the measured gamma-ray energy flux, 
an 100% upper limit for the gamma-ray efficiency and different 
values for /n, including best-fit values obtained from radio and 
gamma-ray light curve modeling (see Section 13.11 and Table [2]l . 
Normally / is assumed to be 10^® gcm^, this would imply a max- 
imum possible distance of about 1.8 kpc for /q = 1. If the distance 
is measured precisely, these limiting curves can be re-interpreted 
as lower limits on the moment of inertia. At the maximum possi- 
ble distance dmax (to be refined as the measurement of the proper 
motion improves) the whole observed Pots would be due to kine- 
matic effects, implying Pcorr = 0. Therefore, for the energy loss 
of the pulsar E = in^IPcorr/P^ to account for the gamma-ray 
luminosity, / would have to be infinite. 



for the small flux density. Furthermore, both pulsars appear 
to be members of very stable MSP-WD binaries. For this 
reason, they have been included in the pulsar timing ar- 
ray (PTA), which is now being used in a collective effort to 
detec t low-frequency gravitational waves (e.g. Hobbs ct al. 
l20ld ). 

Given this high timing precision, continued Arecibo 
timing might provide a precise distance measurement which, 
given the observed gamma-ray energy flux, would provide a 
lower limit for the moment of inertia of this MSP. This dis- 
tance is likely to be significantly smaller than the estimate 
provided by the NE2001 electron model of the Galaxy, in 
which case the parallax should be easier to measure but the 
lower limit on / would not be constraining. If the distance is 
comparable to the DM prediction or larger, then the paral- 
lax will be more difficult to measure, but a low upper limit 
for the parallax would allow us to derive a high lower limit 
for the moment of inertia. This might constrain the EoS of 
neutron matter at densities above that of t he atomic nu- 
cleus - some predict I < 1.7 x 10*^ g cm'^ jWorlev et al.l 
I2OO8I ): measuring a larger lower limit for I would exclude 
such EoSs. 

Continued timing of PSR J2043+1711 will also improve 
the measurement of the Shapiro delay, providing precise esti- 
mates of the masses of the components of the system. Com- 
bining this precise mass with a lower limit on the moment 
of inertia could provide a stringent constraint of the EoS. 
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APPENDIX A: RADIO SEARCH 
OBSERVATIONS WITH THE NANQAY RADIO 
TELESCOPE 

To avoid redundant observations with other instruments, 
and to inform the community of existing data samples. Table 
lAll lists all pulsar searches of Fermi error ellipses made with 
the NRT to date. 

For each observation. Table lAll lists the observation 
time Tabs, the sky temperature in the corresponding direc- 
tion Tsky, and the minimum detectable flux density Smin- 
The sky temperature T^ky corresponds to the contribution 
from the Galactic synchro tron component, cal culated by 
scaling the 408 MHz map of lHaslam et all (j 19821 ) to the ob- 
serving frequency of 1.4 GHz with a spectral index of —2.6. 
The quantity Smin was es timated using the modifl ed ra- 
diometer equation (see e.g. iLorimer fc Kramejiioosi ). with 
G = 1.4 K Jy-\ np = 2, /3 = 1.05, Tsys = Tr,c + T^ky 
with Tree = 35 K (note that this receiver temperature in- 
cludes the 2.7 K temperature from the cosmic microwave 
background). A/ = 128 MHz, and assuming {S/N)jnin = 5, 
and W — 0.1 X P. Also listed in the Table are telescope 
pointing directions, and the offsets from the corresponding 
2FGL sources. 
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For each source, Table lAll gives the semi-major axis 
of the 95% confidence ellipse (S95), the curvature signif- 
icance ("Signif_Curve"), the variability index ("Variabil- 
ity Judex"), and the name of the identified or likely associ- 
ated source ("Assoc."), if any. Details on the determination 
of these parame ters and on the sou rce association procedure 
can be found in lAbdo et all l|2012l ). 

Of the sources observed with the NRT and for which 
no radio pulsations have been observed to date, the follow- 
ing have curvature significances above 4(t and variability in- 
dices smaller than 41.6, making them good pulsar candi- 
dates: 2FGL J0224.0-h6204, J0734.6-1558, J1120.0-2204, 
J1311. 7-3429, J1625. 2-0020, and J2339. 6-0532. With the 
exception of 2FGL J0734.6— 1558, in which a gamma-ray 
pulsar has been d iscovered through bl ind searches of the 
Fermi LAT data (|Saz ParkinsonI uOllI ). these sources re- 
main unassociated. Multi-wavelength studies might help de- 
termine their natures. Optical and X-ray observations of 
2FGL J2339.6-0532 showed that it is likely to be pow- 
ered by a millisecond pulsar in a bl ack-widow system 
l|Romani fc Shawll201ll : fKong et al.ll2012l ). 
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Tabic Al: Fermi LAT sources searched for radio pulsars at the Nangay Radio 
Telescope. Definitions of parameters and details on the calculation of the sky 
temperature T^i^y and the minimum detectable flux density Smin for each 
observation are given in the text. Note that Right ascensions and declinations 
refer to the telescope pointing directions, and not necessarily to the locations 
of the FcTrrri LAT soiircc^s. 



2FGL Name 


IFGL Name 


I 

(deg) 


b 

(deg) 


696 
(arcmin) 


Date 


R.A. 
(J2000.0) 


Dccl. 
(J2000.0) 


Offset 
(arcmin) 


^obs 
(min) 


^sky 
(K) 


(MJy) 


Signif_Curve 


Variability -Index 


Assoc. 


J0031. 0+0724 


J0030.7+0724 


114.10 


-66.11 


7.2 


2010-01-08 


00" 30" 


45= 


+07° 22 


34" 


5.7 


58.5 


0.9 


47 


0.8 


20.9 




J0102.7+B827 


J0102. 8+6827 


124.42 


-4.38 


3.6 


2009-12-17 


0l''02"' 


50= 


+58°27 


18" 


0.8 


66.0 


2.0 


46 


3.1 


119.3 


TXS 0059+681 












2009-12-23 


0l''02"' 


50= 


+ 58° 27 


18" 


0.8 


66.0 


2.0 


46 


















2010-01-09 


0l''02" 


50= 


+ 58° 27 


18" 


0.8 


66.0 


2.0 


46 


















2010-01-10 


0l''02"' 


50= 


+ 58° 27 


18" 


0.8 


78.0 


2.0 


42 








J0118. 8-2142 


J0118.7-2137 


173.46 


-81.73 


2.9 


2010-04-18 


Ol^lS" 


43= 


-21°37 


39" 


5.4 


57.2 


0.9 


48 


4.1 


257.2 


PKS 0116-219 












2010-05-14 


0l"l8°' 


43= 


-21°37 


39" 


5.4 


54.6 


0.9 


49 








J0124.5-0621 


J0124. 6-0616 


145.21 


-67.79 


12.1 


2011-03-30 


01*' 24" 


35" 


-06°21 


61" 


0.2 


29.7 


0.9 


66 


2.7 


28.8 


PMN J0124-0624 












2011-04-01 


0l''24" 


35° 


-06°21 


61" 


0.2 


39.7 


0.9 


68 


















2011-04-30 


01*' 24" 


35° 


-06°21 


61" 


0.2 


49.7 


0.9 


61 








J0127.2+0324 


J0127.0+0322 


140.12 


-68.26 


5.8 


2011-03-26 


0l''27" 


17° 


+03°26 


11" 


0.8 


61.0 


0.9 


61 


0.1 


31.9 


NVSS J012713+032259 












2011-03-29 


0l''27" 


17° 


+03° 26 


11" 


0.8 


61.0 


0.9 


61 


















2011-05-07 


0l''27" 


17° 


+03°25 


11" 


0.8 


62.0 


0.9 


60 








J0131. 1+6121 


J0131. 2+6121 


127.67 


-1.15 


2.5 


2009-12-22 


Ol^Sl" 


17= 


+61°21 


42" 


1.1 


66.0 


2.4 


46 


2.4 


58.1 


IRXS J013106.4+612035 












2009-12-30 


0l''31" 


17= 


+ 61° 21 


42" 


1.1 


58.6 


2.4 


49 


















2009-12-31 


0l"31" 


17= 


+ 61°21 


42" 


1.1 


34.7 


2.4 


64 


















2010-09-25 


()l''31" 


17= 


+ 61°21 


42" 


1.1 


53.0 


2.4 


52 


















2010-12-02 


0l''3l" 


08= 


+ 61° 20 


45" 


0.6 


31.0 


2.4 


68 








J0137.7+5811 


J0137.8+6814 


129.02 


-4.12 


6.0 


2010-01-12 




46= 


+ 58° 14 


07" 


2.2 


34.7 


2.0 


63 


0.7 


15.3 


IRXS J013748.0+581422 












2010-01-13 


0l''37" 


46= 


+58° 14 


07" 


2.2 


66.0 


2.0 


46 


















2010-01-20 


0l''37" 


46= 


+58° 14 


07" 


2.2 


32.5 


2.0 


66 


















2010-12-01 


0l'"38" 


04= 


+58° 13 


12" 


2.7 


23.5 


2.0 


77 


















2010-12-14 


0l''38" 


04° 


+58° 13 


12" 


2.7 


39.7 


2.0 


69 








J0224.0+6204 


J0224.0+6201C 


133.55 


1.13 


3.4 


2011-04-08 


02*' 24" 


28° 


+62°02 


59" 


3.1 


47.4 


3.0 


66 


7.2 


20.1 




J0260. 7+5631 


J0251. 5+6634 


138.87 


-2.61 


7.1 


2009-12-22 


02'' 50" 


47° 


+56° 31 


50" 


0.3 


60.0 


2.1 


48 


0.6 


28.0 


NVSS J026047+562935 












2009-12-30 


02'>50" 


47° 


+56°31 


50" 


0.3 


66.1 


2.1 


60 


















2009-12-31 


02'' 50" 


47= 


+56°31 


50" 


0.3 


32.6 


2.1 


66 








J0334.3-3728 


J0334. 4-3727 


240.22 


-64.36 


3.1 


2009-11-16 


03'' 34" 


32= 


-37° 27 


43" 


2.3 


67.6 


0.6 


47 


2.7 


90.3 


PMN J0334-3726 












2010-01-14 


03'' 34" 


32= 


-37° 27 


43" 


2.3 


60.0 


0.6 


46 








J0434.1-2014 


J0434. 1-2018 


217.85 


-38.94 


10.2 


2009-12-15 


04'' 34" 


13= 


-20° 18 


lO" 


4.0 


52.2 


0.7 


50 


1.1 


42.4 


TXS 0431-203 












2010-04-14 


04'' 34" 


13= 


-20° 18 


10" 


4.0 


19.7 


0.7 


81 


















2010-04-23 


04'' 34" 


13= 


-20° 18 


10" 


4.0 


13.5 


0.7 


98 


















2010-04-24 


04'' 34" 


13= 


-20° 18 


10" 


4.0 


13.5 


0.7 


98 


















2011-04-05 


04'' 34" 


13= 


-20° 18 


lO" 


4.0 


44.7 


0.7 


54 








J0523. 3-2530 


J0523. 5-2529 


228.23 


-29.84 


4.0 


2009-11-08 


05'' 23" 


32= 


-25°30 


26" 


3.0 


32.2 


0.6 


63 


2.1 


22.7 














2009-11-17 


06'' 23" 


32° 


-25°30 


26" 


3.0 


24.7 


0.6 


72 








J0714. 0+1933 


J0714.0+1936 


197.68 


13.61 


2.6 


2009-11-06 


07'' 14" 


02° 


+19°34'39" 


0.8 


38.6 


1.2 


69 


5.1 


313.4 


MG2 J071364+1934 


J0734.6-1558 


J0734.7-1667 


232.04 


2.01 


3.4 


2010-09-07 


07'' 34" 


45° 


-15°69 


18" 


1.4 


66.0 


1.4 


46 


5.9 


24.5 


LAT PSR J0734-1669 


J0928. 8-3530 


J0929. 0-3631 


263.02 


11.23 


7.7 


2010-02-08 


09'' 28" 


56° 


-35°32 


16" 


2.1 


23.6 


1.4 


76 


2.8 


23.2 




J1120.0-2204 


J1119. 9-2206 


276.50 


36.06 


3.8 


2010-01-14 


ll''19" 


54° 


-22°06 


06" 


1.5 


44.7 


0.6 


64 


5.5 


25.9 














2010-01-17 


ll''19" 


54° 


-22°06 


06" 


1.5 


36.0 


0.6 


60 


















2010-07-18 




54= 


-22°05 


06" 


1.5 


60.0 


0.6 


46 


















2011-03-29 


ll''19" 


58= 


-22° 05 


26" 


0.8 


56.0 


0.6 


48 


















2011-03-30 


ll''19"' 


58= 


-22° 05 


26" 


0.8 


57.2 


0.6 


47 


















2011-04-22 


ll''19" 


58= 


-22°05 


26" 


0.8 


47.2 


0.6 


52 








J1120. 4+0710 


J1120.4+0710 


251.53 


60.61 


5.3 


2009-11-27 


11 ''20'" 


30= 


+ 07° 10 


22" 


0.5 


27.2 


0.7 


69 


2.8 


26.6 


MGl J112039+0704 












2009-12-03 


11 ''20" 


30= 


+ 07° 10 


22" 


0.5 


66.0 


0.7 


44 


















2010-07-13 


ll''20" 


30= 


+07° 10 


22" 


0.5 


66.0 


0.7 


44 


















2010-07-14 


ll'"20" 


30= 


+07° 10 


22" 


0.5 


66.0 


0.7 


44 








J1129. 5+3758 


J1129. 3+3757 


175.54 


69.69 


8.7 


2009-11-19 


ll''29" 


23= 


+37° 58 


19" 


1.6 


55.4 


0.8 


49 


0.4 


28.3 
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2PGL Name 


IFGL Name 


I 

(deg) 


b 

(deg) 


"95 
(arcmin) 


Date 


R.A. 

(J2000.0) 


Dccl. 
(J2000.0) 


CMsct 
(arcmin) 


Tabs 
(min) 


Tsky 

(K) 




Signif_Curve 


Variability_Index 


Assoc. 












2010-01-11 


ll''29° 


23= 


+37° 58 


19" 


1 


6 


32 


9 


0.8 


63 


















2010-05-10 


ll''29" 


23" 


+37°68 


19" 


1 


6 


29 


7 


0.8 


66 


















2010-05-12 


ll''29" 


23= 


+37° 68 


19" 


1 


6 


26 


1 


0.8 


71 


















2010-07-21 


ll''29" 


23= 


+37°68 


19" 


1 


6 


39 


7 


0.8 


67 


















2010-12-01 


ll''29" 


56= 


+37°69 


05" 


5 





37 


2 


0.8 


69 


















2010-12-12 


ll''29" 


56= 


+37°69 


05" 


5 





66 





0.8 


44 








J1311. 7-3429 


J1311. 7-3429 


307.69 


28.20 


2.0 


2010-03-30 




47= 


-34°29 


39" 





4 


60 





1.5 


48 


6.3 


19.1 














2010-05-17 


13'' 11" 


47= 


-34° 29 


39" 





4 


37 


2 


1.5 


60 


















2011-05-27 




47= 


-34° 29 


39" 





4 


37 


6 


1.5 


60 








J1340.5-0412 


J1340. 5-0413 


325.51 


66.60 


8.7 


2010-03-28 


13''40" 


33= 


-04° 13 


56" 


1 


5 


60 





0.9 


47 


0.3 


25.0 














2010-04-13 


13'Mo" 


33= 


-04° 13 


56" 


1 


5 


60 





0.9 


47 


















2010-07-21 


13'M0" 


33= 


-04° 13 


5B" 


1 


5 


32 


2 


0.9 


64 


















2010-09-07 


13'M0" 


33= 


-04° 13 


5B" 


1 


5 


41 





0.9 


57 


















2()l()-09-()9 


13'M0° 


33= 


-04° 13 


56" 


1 


5 


41 





0.9 


57 


















2()l()-12-()2 


IS'MO" 


33= 


-04° 13 


56" 


1 


5 


53 


5 


0.9 


50 


















2010-12-04 


13'' 40° 


33= 


-04° 13 


56" 


1 


5 


31 





0.9 


65 


















2010-12-08 


13'' 40° 


33= 


-04° 13 


66" 


1 


5 


27 


2 


0.9 


69 


















2011-02-11 


13*' 40° 


33= 


-04° 13 


66" 


1 


5 


62 


2 


0.9 


60 


















2011-03-22 


13*' 40° 


33= 


-04° 13 


66" 


1 


5 


61 





0.9 


61 


















2011-05-09 


13*' 40° 


33= 


-04° 13 


66" 


1 


5 


44 


7 


0.9 


64 


















2011-05-13 


13*' 40° 


33= 


-04° 13 


56" 


1 


5 


46 





0.9 


63 








J1407.4-2948 


J1407.9-2928 


322.01 


30.22 


21.1 


2011-03-27 


14*' 07° 


25= 


-29°48 


59" 





1 


68 


6 


1.5 


48 


3.2 


17.6 














2011-03-31 


14''07° 


25= 


-29°48 


59" 





1 


55 





1.5 


50 


















2()ll-()4-()2 


14'' 07" 


25= 


-29°48 


59" 





1 


22 


2 


1.5 


78 


















2011-04-27 


14'' 07° 


25= 


-29°48 


59" 





1 


44 


7 


1.5 


55 


















2011-05-07 


14'' 07° 


25= 


-29°48 


59" 





1 


43 


5 


1.5 


56 


















2011-05-11 


14'' 07° 


25= 


-29°48 


59" 





1 


66 





1.5 


45 








J1625. 2-0020 


J1625. 3-0019 


13.92 


31.83 


3.4 


2009-11-11 


16'' 25° 


18= 


-00° 19 


15" 


1 


5 


54 


7 


1.4 


50 


8.3 


24.6 














2009-11-19 


16'' 25° 


18= 


-00° 19 


15" 


1 


5 


46 





1.4 


54 


















2010-01-12 


16''25° 


18= 


-00° 19 


15" 


1 


5 


60 





1.4 


47 








J1642. 9+3949 


J1642. 5+3947 


63.48 


40.95 


3.2 


2009-11-04 


16*' 42° 


32= 


+39° 46 


15" 


6 


3 


66 





0.9 


46 


1.3 


111.7 


3C 345 












2009-11-20 


16'' 42° 


32= 


+39° 46 


15" 


6 


3 


66 





0.9 


46 








J2001. 1+4352 


J2001. 1+4361 


79.06 


7.12 


1.2 


2009-11-06 


20'' 01° 


11= 


+43° 62 


07" 





7 


24 


3 


2.9 


78 


2.1 


118.0 


MAGIC J2001+436 












2009-11-09 


20'' 01° 


11= 


+43°62'07" 





7 


61 





2.9 


64 


















2010-01-10 


20''01° 


11= 


+43° 62 


07" 





7 


67 


4 


2.9 


61 








J2339. 6-0532 


J2339. 7-0631 


81.36 


-62.47 


2.5 


2009-11-21 


23'' 39° 


44= 


-05°31 


13" 


1 


9 


11 


9 


0.8 


105 


5.7 


15.7 














2009-11-21 


23'' 39° 


44= 


-05°31 


13" 


1 


9 


21 


8 


0.8 


77 


















2009-11-26 


23'' 39° 


44= 


-05°31 


13" 


1 


9 


46 





0.8 


63 


















2010-01-12 


23'' 39° 


44= 


-05°31 


13" 


1 


9 


66 





0.8 


48 









